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' Abstract 

p i' We use an analogy between the Yang-Mills theory Hamiltonian and the matrix model 

<D '. 

^ . description of the closed bosonic membrane theory to calculate the spectrum of glue- 



balls in the large Nc limit. Some features of the Yang-Mills theory vacuum, such as 
the screening of the topological charge and vacuum topological susceptibility are dis- 
cussed. We show that the topological susceptibility has different properties depending 
on whether it is calculated in the weak coupling or strong coupling regimes of the 
theory. A mechanism of the formation of the pseudoscalar glueball state within pure 
Yang-Mills theory is proposed and studied. 
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Introduction 



Gluodynamics [Q, being the asymptotically free theory of colored massless 
vector particles, is believed to underline the dynamics of strong interactions. Because 
of asymptotic freedom the theory is well studied at short distances, however, long 
distance phenomena deserve to be understood much better. 

The theory predicts glueballs [Q, the nonperturbative bound states composed of 
pure glue ||^. That prediction was confirmed some time ago by "observing" glue- 
balls in lattice QCD simulations [0, [0. In addition to that, there are experimental 
signatures of resonances which strongly resemble properties of glueballs (for a recent 
analysis of these issues see ref . [|] ) . 

Studying glueballs one might hope to learn more about the complicated ground 
state structure of non-Abelian Yang-Mills (YM) theory. The main question one might 
wonder about is the mechanism of the formation of glueball states in YM theory. 
Those states appear to be heavy in comparison with the lightest hadrons and range, 
depending on the spin-parity structure, within the mass interval 1.5 — 2.3 GeV P], 
1^. Thus, the naive picture of the glueball as a system of two massless gluons which 
interchange virtual perturbative gluons does not seem to be appropriate. 

In this work we are looking for qualitative features responsible for the process of 
formation of pseudoscalar glueballs. A possible mechanism will be proposed. As an 
outcome we calculate the spectrum of lightest pseudoscalar glueballs. The results are 
in agreement with predictions of lattice calculations §]. The paper deals with pure 
YM theory, no light fermions are included. A brief discussion of full QCD is given at 
the end of the work. 

Our study relies on the existence of the 9 term in pure YM theory. We de- 
fine the topological charge density operator as Q = gg^G^^^G^^, with G^^, being 
the nonabelian gauge field strength tensor and the dual tensor is normalized as 
~ k^^J.l^Q|3^a|3■ Because of instantons 0, the non-Abelian gauge theory pos- 
sesses a complicated vacuum structure |T0|. That is, there is an infinite number of 
degenerate vacuum states labeled by some topological invariant, the winding num- 
ber or topological charge. Instantons, being defined in Euclidean space, provide that 
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quasi-classical tunneling processes happen between the different vacua. Thus, the 
true ground state of the theory is a superposition of the vacua with different topo- 
logical charges. The superposition can be provided in the path integral formulation 
by adding to the action the 9 term /\Sg = 6 J d^xQ{x) [T^. However, such a modifi- 
cation is not a harmless procedure. The 6 term in QCD leads to an induced neutron 
electric dipole moment. Experimental bounds on that quantity restrict the value of 
the 6 parameter to be unnaturally small, less than the billion'th part of the unity, 
and give rise the famous strong CP problem ||lT |. 

The picture outlined above implies that the integral of the topological charge 
density / d'^xQ^x), being the topological charge in Euchdean space, is quantized if 
the instanton boundary conditions are imposed on gauge fields Q. Thus, the whole 
scenario of the superposition of the different vacua by means of the 6 term relies on 
quasiclassical arguments. In general, one expects that the quasiclassical approxima- 
tion is justified in a weak coupling regime only |jl2[. What the 6 term leads to in a 
strong coupling approximation where quasiclassical arguments cease to be valid, is 
not clear a priori. 



It was argued by Witten |T3| that in the confining phase of the theory noninteract- 
ing instanton boundary conditions should not be relevant. The key observation was 
that gauge fields with instanton boundary conditions do not yield the area law for 
the Wilson loop. Thus, configurations with those boundary conditions, as any config- 
urations tending to a pure gauge at infinity, fail to satisfy the confinement criterion 



14|. As a result, in the strong coupling approximation of the confining theory one 



should rather encounter some smeared distributions of interacting topological charges 
as opposed to the noninteracting instanton system with quantized topological charge 
|I3| . This statement finds support in recent models of the YM vacuum. Properties of 
hadrons are correctly described by the model in which instantons and antiinstantons 



are coupled in moleculelike (or even more complicated) entities [0 . This kind of in- 
stanton clustering, indicating on strong correlations between them, was also observed 



* Under those boundary condition one means that the vector potential tends to a pure gauge 
configuration at spatial infinity U^^{x)d^U{x), with U{x) being an element of the SU{Nc) 

gauge group. 
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in some recent lattice QCD studies p6| . 

Interactions between instantons, if sufficiently strong, lead to the screening of 
the topological charge at finite distances [l^, [T^, in analogy with the well-know 
phenomenon in plasma physics. 

A quantitative study of the screening phenomenon from fundamental principles is 
a cumbersome task. However, as we mentioned already, both hadron phenomenology 
and lattice simulations seem to support that picture. Below in Section 1 we present 
some arguments (other than the ones mentioned above) indicating that the screening 
of the topological charge should really exist. We show (Section 2) that the three-form 
composite field, which is Hodge dual to the Chern-Simons current and is known to 
propagate the Coulomb type interaction, should be relevant for the description of the 
screening phenomenon. 

In Section 2 we present a possible mechanism of the formation of a pseudoscalar 
glueball. The main ingredients needed for that mechanism to be realized in a confining 
theory are the screening of the topological susceptibility and the presence of the 9 
angle in the theory. 

Having this part set, we discuss in Section 3 an analogy between the spectrum 
of the YM Hamiltonian and that of a closed bosonic membrane with the topology 
of a sphere. We use that analogy to calculate the spectrum of glueballs. In fact, 
we derive a matching condition between the spectrum of YM theory in the large Nc 
limit and the matrix quantum mechanics formulation of the closed bosonic membrane 
theory. Then, studying the spectrum of a spherical closed membrane and using the 
matching condition we calculate the spectrum of the YM Hamiltonian. That gives the 
prediction for masses of lightest glueballs. We also show that the scenario discussed 
in this work is realized only when the 9 parameter is a macroscopic number, i.e. a 
number of order of the unity or so. More accurate estimates are given in Subsection 
3.3. We briefly discuss how the strong CP violation, being present in pure YM theory, 
might still not be observable in full QCD with light quarks. 

Discussions in the present paper are based on the results obtained by M. Liischer 



[111, H, H, by J. Goldstone and J. Hoppe (see ref. Q), and by B. de Wit, J 



Hoppe and H. Nicolai [51|. Where it is possible we present below brief summaries of 
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those results. 



1. Topological Susceptibility 

In this section we study properties of the correlator of the vacuum topological 
susceptibility. We work in Euclidean space-time assuming that the theory is defined 
in a compact Euclidean four- volume V = x r, with / being the linear size of the 
volume and r stands for Euclidean time. The correlator of the vacuum topological 
susceptibility can be written as follows: 



xiV)^ / {0\TQ{x)Qmo)d\ (1) 

V 



where Q is the topological charge density operator defined in the previous section^]. 

The function x(^) is determined by the nonperturbative sector of Yang-Mills 
theory. Calculation of x(^) i^i that respect is a matter of modeling of the vacuum 
structure of Yang-Mills theory. This in its turn is a complicated task. Below we show 
however, that one can still study qualitative features of the volume dependence of the 
correlator of the vacuum topological susceptibility for small and for large values of 
the volumeQ Let us first define, following |T^, what one could call a small volume or 



large volume limit. 

There is a dynamically generated mass scale in YM theory, Aym- The reciprocal 
quantity of Aym sets the characteristic correlation length for the model. Let us denote 
that length by C = AyM- Restricting for simplicity the YM beta function to the next- 
to-leading order approximation, the expression for C can be written as follows:[| 

C = fi ^(aj^'o exp ' ' 



where /x is the renormalization scale, = asifi"^ / Ay^) is the scale-dependent strong 
coupling constant, and Pq and Pi are the first two scheme independent coefficients of 



^ The correlator (1) is in general a divergent Green's function. These divergences can be removed 

by means of standard procedures which are discussed in the Appendix. 

*We vary I keeping r fixed. Thus, we actually study the three- volume dependence of x(V^). 

^ Since we are interested in the magnitude of this quantity the explicit 9 dependence in this 

expression is dropped. 
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the beta function, Pq = llNc/3, Pi = SAN^/S. The expression for ( is renormahzation 
group invariant in the corresponding order. 

Let us now introduce the following two limits. One can define the value of the 
volume element V to be small if the correlation length ( is much larger than V^^^, 
i.e., C >> V^^'^ |]19[. The large volume limit in that case would refer to a volume 
element satisfying the condition C << V^^^. 

One can show that the two limits defined above correspond respectively to the 
weak coupling and strong coupling regimes of the theory. In order to see this let 
us keep the product of the renormahzation scale fi and the value of V^^^ fixed, say, 
Vfi'^ = 1. This condition sets the scale /i as an infrared cutoff. Then, the expression 
for the correlation length takes the form 

[as) exp 



Thus, the small volume approximation given by the condition C >> V^^^ corresponds 
to the weak coupling regime, i.e., << 1. For instance, if one sets ^^7^74 ~ 10 >> 1, 
the corresponding value of the coupling constant is as — 0.2|. 

Let us now turn to the discussion of the large volume limit defined by the condition 
( « V^l'^. This limit is equivalent to the small pi, approximation. However, for small 
values of the running coupling constant as is a big number. Hence, the large volume 
limit corresponds to the strong coupling regime of the theory. Regretfully, we can not 
estimate (as we did in the case of the small coupling constant) how large the coupling 
constant should actually be. The approximation which is used defining C, breaks down 
for large values of the coupling. Though one could present the definition of C for any 



orders of perturbation theory (see, for example, pO|), that definition would contain 

the exact form of the beta function P{a^ which is known only perturbatively. So, the 

all order formula also becomes inappropriate for practical calculations in the strong 

coupling approximation. 

Below we show, at least qualitatively, that the topological susceptibility as a 

function of V has different behavior depending whether it is calculated in the weak or 

§If one approximates Aym — (100 — 200) MeV, then the smaU volume hmit refers to V^l^ << 
fl - 2) fm. 
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strong coupling regimes. In the weak coupling phase it is an increasing function of the 
argument, and on the contrary, in the strong coupling regime the function decreases 
with the argument monotonically. 



1.1. The Weak Coupling Approximation 

Let us start with the small volume or weak coupling approximation. The non- 
Abelian gauge theory provides a good description of physics in that domain. Excita- 
tions with a zero topological charge do not contribute to the value of x(^) defined 
in eq. (1). Only nontrivial topological configurations of gauge fields are to be taken 
into account. In the weak coupling regime the YM vacuum can be approximated by 



noninteracting, well-separated instantons |2l|. In that approximation instantons can 



be treated as point-like objects. The expression for the topological charge density 
can be written as follows: 

i 

where denotes the topological charge for a configuration localized at the point Xj. 
Assuming that instantons do not interact with one another we derive 

1 +00 

(0|TQ(x)g(0)|0) = - E ^^"\^), (2) 

q= — oo 

where the index i in the definition of the topological charge was omitted. The quan- 
tity Pq denotes the probability for a nonabelian gauge field configuration to have a 
topological charge equal to q. These probabilities are exponentially suppressed for 
nonzero q and one expects that the infinite series in eq. (2) converges^]. 
Substituting eq. (2) into eq. (1) one derives 

1 +00 

x{v) = y y: ^"p.iv)- (3) 

In the approximation we set above the following relation is valid: 

The partition function Z{9) in that case can be approximated as Z{9) = J^q Pq Gxp{i9q). 



where Pi{V) denotes the one instanton contribution. Substituting this relation into 
eq. (3) and performing the summation of the infinite series one gets the following 
expression for the topological susceptibility: 

X{V) = ^^^(^^ ^ + ■ (4) 

Thus, the small volume behavior of the function x(^) is approximately defined by 
eq. (4). The expression for P\{y^ can be calculated in the one loop approximation 
using the well known results of ref. 

P,(l-) = const xexp ['^J^^^TJ^ ^ 
where the following expression for the strong coupling constant is supposed to be used 

The result for -Pi(l^) is 



/?oln(Vi/2A^ 
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P\{y) = const X 1^1 X logarithms, 

where the logarithmic corrections appear in the next-to-leading approximation. Hence, 
as y — the ratio ^-^y-^ also tends to zero. As a consequence, in the small volume 
limit limy^ox(^) ~^ 0. Moreover, based on the relations given above one concludes 
that for small volume elements the quantity x(^) is a monotonically increasing func- 
tion of the argument V. This property should hold as the condition C >> V^^^ is 
satisfied. 

Suppose now that the quantity V^^^ becomes comparable in magnitude with ( so 
that the weak coupling approximation breaks down. As a result, the pointlike nonin- 
teracting instanton approximation ceases to be valid. Interactions between instantons 



start to play a crucial role providing the screening of the topological charge [T^ . 

Let us assume for a moment that one neglects instanton interactions even for large 
values of the volume and let us study what happens in this unrealistic case keeping in 
mind that the interaction effects are going to be included later. Doing so one is dealing 
with an ideal gas of instantons placed in a large volume. The approximate calculation 
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of the partition function with noninteracting instantons in the thermodynamic hmit 
yields the following Gaussian distribution function for PqiV) 



V2nVd V 21^^/ 

where ci is a not yet defined positive constant. We substitute this expression into eq. 
(3) and perform the summation of the infinite series in the large volume limit. The 
final expression can be found using the following relation 



1 +00 



■wq^/b 

Vbq-1 



where b is an arbitrary positive number . As a result one gets 



lim xiV) = d. (5) 

Let us summarize briefly what we learned about the volume dependence of the 
topological susceptibility x(V"). In the zero volume approximation the topological 
susceptibility was zero. Increasing the value of V , so that the weak coupling approx- 
imation still holds, the function x(^) increases monotonically. If one goes further 
and neglects the interaction between instantons even in the large volume (strong cou- 
pling) approximation, one finds that the function xiV) reaches its asymptotic valueQ 
denoted above by xiY >> C^) = d. However, as we stressed earlier, interactions be- 
tween instantons play a crucial role in the strong coupling approximation. In the next 
subsection we show that the topological susceptibility becomes a decreasing function 
of the argument for large values of V when the effects of finite distance correlations 
between topological charges are taken into account. 

1.2. The Strong Coupling Approximation 

Let us consider the large volume or strong coupling limit. In that limit the theory 
is in a confining phase. Instantons are interacting strongly^. Those interactions 



^ These properties were originally studied in ref. p9[ | considering YM theory on a four-sphere S^. 
^li is not even clear whether it makes sense to talk about a configuration with a definite topological 

charge in this case |13[. 
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become responsible for formation of spin zero glueball states ||2^. A description in 
terms of colored variables is not a good approximation anymore. The theory, however, 
can be defined by means of the low-energy effective action containing colorless degrees 
of freedom. The explicit form of that effective action for pure YM theory is not known. 
In general, the action can be written as 



Seff = J d'xC{Gn,VGn,V^Gn,...), 



where G^s stand for glueball fields. 

Below we study properties of the correlator of the vacuum topological suscepti- 
bility in the effective theory. We denote this quantity by XcfT(^)- The correlator in 
eq. (1) is saturated by the set of intermediate glueball states 

{0\TQ{x)Q{0)\0) = d 5^-^\x) + Y.{^\Q\n){n\Q^Dp{mn\x\), (6) 

n 

where m„ is the mass of the n th intermediate physical state and Di?(m„|x|) stands 
for the Euclidean x-space Feynman propagator of a scalar massive particle 

DF{jnn\x\) = Ki{mn\x\), 

with Kilrrinlxl) being the Bessel function of an imaginary argument. 

The parameter d given in eq. (6) is a positive number. It was introduced in the 
preceding section and in the simplest case of a dilute instanton gas approximation 
corresponds to the value of the topological susceptibility in the large volume limit. 
From the point of view of the effective theory we deal with, the parameter d is 
a momentum independent subtraction coefficient in the dispersion relation for Xes 
written in momentum spaceQ. In eq. (6) we implicitly assumed that the volume 
element is sufficiently large so that the YM topological susceptibility occurring as the 
first term on the r.h.s. equals to its asymptotic value d. 

Strictly speaking, there are additional continuum contributions on the r.h.s. of eq. 
(6). They account for possible many-particle intermediate states. Those contributions 
^There is another subtraction term in eq. (6). It is proportional to the second derivative of the 
Dirac delta function. This term, being integrated in eq. (1) gives a vanishing contribution and does 
not appear in the definition of x(^)- ^ detailed discussion is given in the Appendix. 
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are studied in the Appendix. We just mention here that the continuum contributions 
do not affect the physical picture we are going to discuss in this subsection. 

One notices that eq. (2), which includes only noninteracting instanton effects, 
reffects the lack of finite distance correlations between topological charge densities, 
i.e. the r.h.s. of eq. (2) is zero for any nonzero value of x. This would not be the 
case if instanton interactions were taken into account. We also saw that the inser- 
tion of the intermediate glueball states into eq. (2) yields the expression (6) with 
finite distance correlations occurring on its r.h.s. Thus, one argues that the strong 
correlations between instantons, which are responsible for finite distance effects, are 
phenomenologically included in eq. (6) as the intermediate glueball states are taken 
into account. The argument above becomes more sensible if one recalls that inter- 
actions between instantons are responsible for the formation of those intermediate 



glueballs p5 |. 



Let us now define the matrix elements occurring in eq. (6). The operator of the 
topological charge density Q is an antihermitian operator in Euclidean space. Taking 
this into account one introduces the following parametrization for the matrix elements 

{0\Q\n) = -lUml (0|Q|n)(n|Q|0) = -^X, 

where /„ can be thought of as a decay constant of the corresponding n'th glueball 
state (in analogy with the pion decay constant /tt). If one substitutes these definitions 
back into eq. (6) the following expression emerges 

(o|rg(a;)g(o)|o) = d s^'\x) -j:fnK^M^n\x\). (7) 

Having this relation established let us study what happens with the correlator of the 
topological susceptibility (eq. (1)). Substituting eq. (7) into eq. (1) we find 

Xefr(V^)=rf-E/nX^n(^), (8) 
n 

where 

Qn{V) =ml DF{'mn\x\)d'^x. 

J \/ 
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The function QniV) determines the volume dependence of the topological susceptibil- 
ity for large values of V. This function has a simple behavior|. The straightforward 
calculation yields 



Qn{0) = 0, 



Qn{00) = 1. 



In general, Gn{V) is a monotonically increasing function of the argument. Its value 
increases rapidly from zero at = to almost its asymptotic value at some finite V. 
Then, increasing very slowly, the function approaches the unity as V — oo. Relying 
on these properties one derives 



Thus, we see that the topological susceptibility gets additional positive subtraction 



Finally, using eq. (8) and the monotonicity of the function Qn{V) one concludes 

that the topological susceptibility decreases from its value defined at relatively small 

volumes to its value reached in the large volume limit. Physically this can be thought 

of as following. Suppose we set a sequence of subvolumes enclosing some topological 

charge distribution Vi < V2 < V3 < ... < oo. The result of our discussion is that 

§ For a spherically symmetric volume element with the radius R one can calculate that Gn{V) — 
l-mlR^K2{mnR)/2. 

^^In general, the quantity on the r.h.s. of eq. (9) is a nonzero number. However, it was argued 
in ref. fv^ that the topological susceptibility might completely be screened in the infinite volume 
limit if instanton interactions are sufficiently strong, i.e. Xcs{oo) would equal to zero in that case. 
This condition would yield a relation between the quantity d and parameters of glueballs. Imposing 
the condition Xoff(oo) — [0 one derives d = J2n fn^'n- This relation is the analog of the Witten- 
Veneziano formula for the 77' meson mass (if one considers full QCD and combines the relation 

derived above with the Witten-Veneziano formula one necessarily needs to take into account the fact 
that the value of d depends on whether it is calculated in pure YM theory or in full QCD). Thus, 
the relation d = fn''^n ^ phenomenological criterion of the validity of the proposal of ref. ||l7|| . 
That relation can be tested in lattice QCD studies by measuring d in noninteracting instanton gas 
picture of pure YM theory and also by studying masses and decay constants of the whole tower of 
pseudoscalar glueball states. For the mass and decay constant of the lightest glueball QCD sum rule 
results can also be used [|8), p9t . 




(9) 



n 



terms in the effective theory (the sum on the r.h.s.)^]. 
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Xefr(^i) > Xefr(V2) > XesiVs) > ...Xcfr(oo). Thus, for smaller volumes one gets larger 
values of the topological susceptibility. One should remember, however, that this 
picture emerges in the confining phase of the theory, i.e., when Vi >> and effective 
degrees of freedom are colorless excitations. The interpretation in terms of colored 
gluons does not make sense in that region because of the lack of asymptotic in and 
out states for those excitations. 

Recalling that the behavior of the function Qn{V) is governed by exponents of the 
type e^~"^"-^\ one concludes that the effective size at which Xefr(^) gets substantially 
suppressed is defined by the Compton wavelength of the lightest 0^+ glueball state 
(eq. (8)). In accordance with lattice calculations the lightest pseudoscalar glueball of 
pure YM theory is expected to have mass approximately equal to 2.3 GeV Q. Thus, 
the effective suppression length scale is L = l/mco — 0.09 fm. 

We complete this subsection by listing the main qualitative conclusions of the 
discussion presented above. 

(i) In the domain of asymptotic freedom, where YM theory is defined most accu- 
rately, the topological susceptibility is an increasing function of the argument. 

(ii) In contrast, in the phase where composite, colorless excitations are formed the 
topological susceptibility decreases monotonically. 

(iii) The suppression length of the topological susceptibility is defined by the 
inverse mass of the lightest 0""*" glueball state and equals approximately 0.09 fm. 
This length is less than the effective radius of the ^ glueball itself (approximately 
O.T-l.OfmUl, H). 

An underlying nonperturbative mechanism which is responsible for the formation 
of the glueball state, most likely, should also be responsible for the suppression 
of the topological susceptibility and vice versa. However, the dynamical reason un- 
derlying these properties is not captured by the qualitative discussion of this section. 
The question why can this happen will be addressed below. 

Let us notice that the properties listed above are in analogy with what happens in 
the (2+l)-dimensional Polyakov [31] model. In the case of the Polyakov model those 



features can be derived in a rather model-independent way [^, 
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2. The 6 Angle and Formation of Glueballs 



In this section we study how gluebaUs can be formed in the vacuum of YM theory. 
In order to address this question let us first recall how quark containing hadrons are 
formed in QCD |Q, |3^. It was found in ref. |2l| that nonperturbative fluctuations 



lower the value of the vacuum energy density in QCD: If the ground state energy 
density for a perturbative vacuum was zero, then instantons lower it yielding a nega- 
tive value . When colored quarks are submerged in that vacuum the QCD ground 



state responds to the insertion of the quarks by suppressing the instanton density in a 
small domain around the quarks [^. In other words, quarks, being submerged in the 
YM vacuum, yield a positive energy density which in the domain around the quarks 
partially compensates the existed negative ground state energy density. The size of 
that domain is determined by the dynamics of nonperturbative QCD [Q. Hence, if 
one takes the value of the vacuum energy density inside the quark containing domain 
and subtracts the value of the vacuum energy density outside the domain one would 
be left with a positive energy density excess in the interior of the domain. 

Having a positive energy excess inside of some region means that there should be 
an inward pressure acting on each small volume element of that domain. In other 
words, the outside region with the negative energy density produces a pressure on 
the quark containing domain tending to squeeze its volume down to zero. The quark 
conflnement emerges in this picture as an effect of the complicated structure of the 
QCD ground state. This serves as a derivation of the bag model for hadrons , 
from fundamental principles of nonperturbative QCD p3[ . 

In this section we show that the same phenomenon might occur in pure YM 
theory. The crucial difference from the previous case is the existence of a purely 
gluonic domain with a positive energy density excess. That positive energy density 
can be provided by the 6 term. We show below that the positive energy density in the 
interior of the domain is proportional to the value of 6^ multiplied by the value of the 
topological susceptibility. Since the topological susceptibility is screened outside of 
some region, this naturally yields a compact region of space with a positive vacuum 
energy density excess inside. We show that this domain can hadronize forming a YM 

14 



glueball state. 

Let us start with the action for Yang-Mills theory with the CP violating 6 term. 
In this section we work in Minkowski space-time. 

We decompose the action S into the usual CP even S^^\ and CP odd S^^\ parts. 

5 = 5(+) + 5(-), where 

s^^^ = J d'x l^-^^Gi,Gi)i , s(-) = ejd'x . (10) 

The total energy density of this system S is the sum of the energy densities of the 
CP even part £^+\ and CP odd part £^-\ i.e. £ = + £'^-\ 

Let us consider the CP even part of the action. As we mentioned above, nonper- 
turbative contributions yield a negative vacuum energy density |]2T]]. The total energy 
density of the CP even part is a sum of the negative ground state energy density £^^ 
and the energy density of perturbations about that ground state f^pg^ 

£W=fW+fS. (11) 

Suppose we start with no perturbations being excited, i.e. put £pcrt = 0- Then, 
S^^^ = ^itJ = i(0|eW|0) = ^{0\d^Gl^\0) ^ -(0.250 GeV)^ [0. Here, 
stands for the anomalous trace of the energy-momentum tensor corresponding to S^^^ 
and perturbative contributions to the gluon condensate are subtracted. 

Let us now address the question what is the contribution of the CP odd part of 
the action to the total energy density of the whole system given in eq. (10). 

It is convenient to introduce a new variable by rewriting the expression for the 
topological charge density^ Q in terms of a four-form field F^^""^ 

where the four-form field F^^"^ is the field strength for the three-form potential C^j^ua 

* Though in Minkowski space-time Q does not have the meaning of the topological charge density 
and, moreover, differs from Euclidean definition of the topological charge by i, we formally keep that 
name and letter for simplicity. 
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The C^ya field itself is defined as a composite operator of colored gluon fields A" 

with fabc being structure constants of the corresponding SU{Nc) gauge group. The 
right-left derivative in this expression acts as AdB = A{dB) — {dA)B. 

The topological charge density can also be expressed through the Chern-Simons 
current as Q = d^K^. Using this expression one can deduce the relation between 
the Chern-Simons current and the three-form potential CyaP) these two quantities are 
Hodge dual to each other K^^ = ■^e'^'^'^^ 0,^0(3 ■ 

Let us rewrite the CP odd part of the action in terms of the three-form potential 
C'i/o/3- For the further convenience the integration over space-time will be restricted 
to a finite, not yet specified domain denoted by A4 

= e [ Qd^x = --, I F.^ap dx^ A dx"" A rfx" A dx^ = -6 [ F, 
Jm 4! Jm Jm 

where the following differential four-form was introduced 

F = ^Ff,^^^ dx^ A dx" A dx"" A dx^. 

In terms of differential forms and an exterior derivative d the equations above formally 
simplify^]. Indeed, F = dC, where C = ^C^ai3 dx" A dx" A dx^^ and the expression 
for S*^^) reads as 

S^-^ = -e I F = -e I C = I C^^B dx'' A dx"" A dx^^. (12) 
Jm JdM 3! JdM 

In the last equation we used the Stokes theorem assuming that the boundary dAi 

enclosing the domain Ad is an orientable smooth surface. Speaking in terms of the 

C^af3 field, the nonzero value of the 6 angle corresponds in Minkowski space to the 

nonzero coupling of the C^ajS field to the boundary manifold dA4. That coupling 

is gauge invariant, although the C^a/s field itself is not a gauge invariant quantity. 

Indeed, if the gauge transformation parameter is denoted by A", then the three-form 

MWe apologize for using the same letter d for an exterior derivative utilized in this Section and 
the quantity d, which has to do with the instanton charge density, defined in the previous Section. 
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field transforms as Cua/B Cvafi + d^h^ap — da^yp — dpKav, where Kap oc A'^dpA"' — 
Al^daA"-. However, it is easy to check that the expression for the field strength F^i^a/j 
is gauge invariant. Since the coupling of the three-form field to the boundary can be 
expressed in terms of the F^ua/s field (as in eq. (12)), then that coupling is also gauge 
invariant and can lead to some physically observable results. The same conclusion 
could be drawn without referring to the field strength. The gauge variation of the last 
expression in eq. (12) is zero for any smooth closed surface which does not enclose 
any field singularities. 

It was noticed some time ago that the Cuap field propagates in the bulk of the 



domain Ai if the topological susceptibility is nonzero in that domain. This becomes 



more evident if one recalls the notion of the Kogut-Susskind (KS) pole in the 
correlator of two Chern-Simons currents. We briefiy present those arguments. 

Consider the correlator of the vacuum topological susceptibility at a nonzero mo- 
mentum. The topological charge density Q is the derivative of the Chern-Simons 
current Q = d^K^. One can substitute this definition back into the expression for 
the correlator of the topological susceptibility. In that way one discovers that x is 
defined as the zero momentum limit of the correlator of two Chern-Simons currents 
multiplied by two momentaj^ 

X = zlimgV / e^«^(0|Tir„(x)ir^(0)|0)ci^x. (13) 

|?->0 J 

The only way for this expression to be nonzero is to claim that the correlator of two 
Chern-Simons currents develops a pole as the momentum goes to zero. This is called 
the Kogut-Susskind pole ||3^ . 



Knowing that the correlator of two Chern-Simons currents has a pole, one can 



^The multiplier —i appears in the definition of tlie topological susceptibility in Minkowski space- 
time. There are some delicate issues regarding the definition of the correlator of the vacuum topolog- 
ical susceptibility. If one defines x ^is a second derivative of the partition function w.r.t. the angle, 
then some contact term appears in that expression [ p6| . Likewise, a special care is needed while 
treating the covariant T product in eqs. (1) and (13) when this last is taken to be the definition of 
X [0- One should add a contact term (given in the Appendix of rcf. [26]) to the r.h.s. of eq. (1) in 
order to define x as a second derivative of the vacuum energy w.r.t. the 9 angle [26,41] (this contact 
term can effectively be included in eq. (6) by redefining the value of the positive constant d). 
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use the relation between the Chern-Simons current and the three-form Cuais field 
and conclude that the C^ai3 field also has a nonzero Coulomb propagator Thus, 
the Ci^ai3 field behaves as a massless collective excitation transferring a long range 



interaction 3S 



Let us summarize briefiy the results of the discussion given above. Following ref. 
[ p8| we established that the three-form field Ci,ai3 propagates in the bulk transferring a 
long-range Coulomb interaction. The exact propagator of this field is of the Coulomb 
type and is proportional to the value of the vacuum topological susceptibility. 

We also saw that the CP odd term in the action of Yang-Mills theory can be 
expressed as a coupling of the three-form composite field Cua/s to the boundary man- 
ifold. Hence, the three-form field Ci^afs being a free field in the bulk actually does 
couple to the boundary surface. 

All the properties mentioned above can be summarized in the following effective 
action for the C^a/s field: 

= -^1^) L - 1 L ^ ' <"> 

The first term in this expression yields the correct Coulomb propagator for the three- 
form C,jap field. The second term is just the usual CP odd 9 term of the initial 
YM action. Vm denotes the three- volume of the domain M.. Notice that higher 
derivative terms are neglected in this action as they are suppressed by momenta of 
massless three-form field. 

Our next step is to study the effective action given in eq. (14) In particular, 
we will calculate the ground state energy of the system using the effective action 
(14). Before we turn to that calculation let us mention that Maxwell's equations for 
a free four-form potential -F^j/a/? yield only a constant solution in (3 + l)-dimensional 
space-time |H3[. The reason is the following. A four-form potential has only one 



^^^One should notice that the action (14) is not an effective action in the Wilsonian sense. It 
is rather related to the generating functional of one-particle-irreducible diagrams of the composite 
field. The effective action in eq. (14) is not to be quantized and loop diagrams of that action are 
not to be taken into account in calculating higher order Green's functions. The analogous effective 
action for the CP even part of the theory was constructed in refs. El[, Ga. 



18 



independent degree of freedom in four-dimensional space-time, let us call it S. Then, 
the four Maxwell's equations written in terms of the S field ensure that this field 
is independent of all four space-time coordinates. Hence, the solution can only be 
a space-time constant. Thus, the F^i,af3 field propagates no dynamical degrees of 
freedom|U]. However, this field can be responsible for the existence of a positive vacuum 
energy density in different models of Quantum Field Theory (see ref. ^4[)- Thus, 
studying classical equations of motion for the F^^q,/? field one can determine the value 
of the ground state energy given by configurations of F^^a/s- We are going to solve 
explicitly classical equations of motion for the effective action (14). Then, the energy 
density associated with those solutions will be calculated. 

Let us start with the equations of motion. Taking the variation of the action (14) 
with respect to the field one gets 

d'^Fnuaaiz) = 6x{Vm) / S^'^^z — x) dx^ A dxa A dxa. 

JdM 



It has been shown in ref. that this kind of equations can be solved exactly in 
four-dimensional space-time. The solution is the sum of a particular solution of the 
inhomogeneous equation and a general solution of the homogeneous equation 

F^ual3{z) = Ox{Vm) / S'^^Hz - X) dXf, A dx^, A dXa AdXp + h Sf^^afS- 

Jm 

The integration constant h, if nonzero, induces an additional CP violation beyond the 
existed 6 angle . Since we would like to avoid to have an extra CP violating term 
we set h = 0. Simplifying the previous equation one finds that the classical solution 
is a nonzero constant tensor density inside of the domain A4 

F^uaP = —dx{^M)^tJ.val3-, (15) 

and is zero, F^jyo/? = 0, outside of M.. 

As a next step let us compute the vacuum energy associated with the solution 
given in eq. (15). The density of the energy-momentum tensor for the CP odd sector 
of the theory can be written as 

p(-) ( p aPr _ 1„ p2 



"The exception is when that field couples to other fields. 
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Using the expression (15) one calculates the corresponding energy density[^ S'^ ^ 

Since the F^j^a/j field does not propagate dynamical degrees of freedom the expression 
above is the total energy density of the system given by the action (14). The crucial 
thing about this energy density is that it is a positive quantity proportional to 6'^ 
multiplied by the value of the topological susceptibility|. 

We learned in the preceding section that the magnitude of the topological sus- 
ceptibility depends on the value of the subvolume in which it is calculated, and also, 
most importantly, it is screened by nonperturbative effects of YM theory outside of 
some finite subvolume element. 

So far we treated the domain Ai as some arbitrary volume. Let us now suppose 
that A4 is the subvolume outside of which the topological susceptibility is screened^. 
Thus, the topological susceptibility is given by the quantity d inside of the volume 

and by the difference d — J2n fn^n outside of the volume. As we mentioned 
above, the difference d — J2n fn^n turns to zero if interactions between instantons are 



sufficiently strong ||T^. However, this is not guaranteed in general. Trying to deal 
with the most general case, we assume here that d — J2n fn^'n is some number not 
necessarily equal to zero. Clearly, our result presented below will also be applicable 
to the case when the topological susceptibility is zero outside of the volume and 
d = J2n fn^n- accordaucc with our previous calculations, the vacuum energy 



**One can use either £ — jO^^ or £ — 9oo- 
■I- One might wonder whether the same result is obtained if one treats 9 not as a constant multiply- 
ing Q in the Lagrangian, but as the phase that the states acquire under a topologically non-trivial 
gauge transformations. These two ways of presenting the 6 dependence are equivalent. Thus, results 
of our discussion should be equivalent in both cases. The key observation is that if 6 is not entering 
the Lagrangian, the arbitrary integration constant in eq. (15) has to be nonzero. It should be chosen 
in a way that would guarantee a proper 9 dependence of the VEV of the topological charge density. 

That would leave the results of our discussion without modifications. 

tin this discussion we assume that the domain has a more or less definite boundary, or in other 

words, that there is a narrow interval where the topological susceptibility drops in magnitude from 

the value d to d~J2n fn^'n- That assumption does not seem to be far from the reality if one recalls 

the behavior of the function Qn{V) (see Section 1). 
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density inside of the domain M. is ^^inside = + \0'^d. Let us now turn to the 
vacuum energy density outside of the volume M.. As we assumed, the topological 
susceptibility is not necessarily zero outside of M.. Hence, the three-form field can 
propagate in that region too and Ff^^a/s 7^ outside of Ai in the general case. As a 
result, there exists some nonzero vacuum energy density outside of Ai. In analogy 
with the previous case one derives the following expression for the vacuum energy 
density £^outside — + — J2n flj^V)- This value is less than the energy density 
inside of the volume M.. Thus, there is a positive energy density excess inside of M.. 
The expression for the corresponding energy excess in the subvolume takes the form 

= \e^^xVM- (16) 

Here Ax = Xinsidc — Xoutsidc = Hn fn^n difference between the topological 

susceptibility defined inside and outside of the subvolume. As we mentioned, eq. (16) 
is also valid in the particular case when the topological susceptibility equals to zero 
outside of the volume M.. In that case Ax = J2n fn'^'i. — d and the energy difference 
in eq. (16) coincides with the energy of the CP odd part inside of the volume M.. 

There are two basic questions to be elucidated here. First of all why would 
any distinguished subvolume exist in the YM vacuum? The reason, as we already 
mentioned above, is the screening of the topological susceptibility. This naturally 
provides finite domains in the YM vacuum with the positive vacuum energy excess 
in accordance with eq. (16). 

The second question is what happens with this finite volume if it is allowed to 
flow freely to a stationary state. The system will tend to minimize the energy given 
in eq. (16). The expansion is not an energetically allowed process. An alternative 
possibility for the system is to squeeze its volume down as much as it is possible. In 
that case the r.h.s. of eq. (16) would be decreasing. In other words, there should 
be an inward pressure acting on the system and tending to squeeze its volume down. 
That pressure is due to the positive difference between the energy densities in the 
interior and exterior of the domain. 

Hence, one concludes that the system will tend to minimize its energy by squeezing 
its volume down, or decreasing Vm in eq. (16). 
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At first glance such a system is unstable and should collapse to a point. However, 
that would be a wrong conclusion. The point is that we did not yet take into account 
perturbations of YM fields which should get excited inside of Vm while the system 
is shrinking its volume down. Those excitations could stabilize the system. The 
energy of those excitations, identified in cq. (11) as S^'^^^, provides an additional 
contribution to the total energy. Anticipating the results of the next section we 
present the expression for the total energy inside of the domain M. 

A£„. = lo^AxV^ + (PO^'t've^numberj ^ ^^^^ 

where the first term is related to the CP odd part of the initial action and the second 
one is the contribution of perturbations of the CP even part. 

The structure of this equation allows one to minimize the quantity A£^totai(KM) 
with respect to Vm and find an optimal value of the three- volume occupied by the 
system. 

We treat this physical system as a model for a pure YM glueball state. 
The spectrum and some properties of that system are studied below. 

3. Yang-Mills vs. Closed Membrane Spectrum 

It was shown in the previous section that the screened topological susceptibility 
leads to a positive energy density excess inside of some finite volume. The system 
tends to minimize that volume. The compression of the volume continues until some 
YM states are exited inside of that domain. Those states have nonzero energy, i.e., 
{'Hym) 7^ 0, where Hym is the Hamiltonian density of YM theory (the CP even 
part) . We mentioned already that the phenomenon described above is related to the 
fact that X is screened. On the other hand, we saw that the effects responsible for 
the screening of the topological charge should also be responsible for the formation 
of the 0^^ glueball state. Hence, it is reasonable to identify those finite volume YM 
excitations with physical ^ glueball states. 

In this section we study the spectrum of the physical YM Hamiltonian in a finite 
volume. Under some approximations, elucidated below, the spectrum of YM theory 
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resembles that of a closed bosonic membrane with the topology of a sphere or torus 
^0| , |51[1 . One can use that analogy to derive the relation between the spectrum of YM 



theory and that of a closed bosonic membrane. Using that relation and calculating 
the spectrum for a closed spherical bosonic membrane we predict masses for two 
lightest pseudoscalar YM glueball states. 

3.1. Studying the spectrum of Yang-Mills theory 

Let us start with the physical Hamiltonian density of YM theory. In order to 
stress the approximations we make we present the brief discussion of the Hamiltonian 
formalism of the theory (for detailed discussions see the textbooks |4^, and Wh )- 



One starts with the Lagrangian density of pure YM theory in Minkowski space- 
time 

The canonically conjugate momentum is defined as the derivative of the Lagrangian 
density w. r. t. the time derivative of the canonical coordinate and is given by 

pa ^ 

9 

The Hamiltonian density takes the form 



9' 1 

£_ pa pa _| 

2 ' ' %2 

2,j = 1,2,3. 



U " pa pa J /~ia /-ia , Aa(T\ r>\a 

/to — —r'i -Tj + ^j— ^(j-jjCTjj + /iQ[Uil^i) , 



This is not the physical Hamiltonian density yet. There are extra degrees of freedom 
in this expression. The existence of those extra variables is related to the gauge 
invariance of the theory. 

The Lagrangian density does not contain time derivatives of the ^0 field. As a 
result, the following primary constraint appears Pq = 0. 

Introduce |^ the so called "total" Hamiltonian Hrit) = J (Px{T-Co + X°'Pq), where 



A"(x) denotes a Lagrange multiplier. Time evolution of a physical quantity is given by 
the Poisson brackets of Ht and the quantity itself. Thus, one needs to set conventions 
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for the Poisson brackets. For any two (bosonic) functionals A and B we use the 
following expression: 

6A 6B 6B 6A \ 



{A, B} = J d^z 



^Sq{z) 6p{z) Sq{z) 6p{z) J ' 

where q and p denote canonical coordinates and momenta respectively. Using this def- 
inition one finds that the conservation of the primary constraint {Pq (x, t), HT{t)} = 0, 
leads one to the secondary constraint in the form of the Gauss's law Df^'Pj' = 0. One 
can also check that the conservation of the Gauss's law is identically satisfied and no 
further constraints are produced at this stage. 

We are going to work in the axial gauge = 0. Requiring the conservation of the 
gauge condition {^3(0;, t), ifr(i)} = 0, one derives the additional secondary constraint 
^2pa _ g^j[a _ Q Finally, the conservation of that constraint leads to the equation 
for determination of the Lagrange multiplier d^X'^lx) + d^diA^lx) — /'^'"^AqSsAq = 0. 
Thus, the whole system of gauge conditions and constraints can be summarized as 



<f ? = PS, $5 = DfP, , 

^ = A^, $^ = g'P^ - d^Al 

The physical Hamiltonian in the axial gauge can be written in terms of the following 
physical variables P^ and A^, where m = 1,2 In general, the straightforward 

procedure implies the elimination of all nonphysical variables by solving (wherever 
it is possible) constraint equations and substituting those expressions back into the 
formula for the Hamiltonian. In most cases the result is a complicated nonlocal 
expression for the Hamiltonian. There is a formally simpler way to follow, however. 
One can solve only some part of the constraint equations keeping the rest of the 
constraints unsolved and allowing some of nonphysical variables to be present in the 
Hamiltonian. Then, the physical system is defined by that Hamiltonian accompanied 
by unsolved constraint equations imposed on the physical states of the Fock space. 
For our purposes we found it convenient to follow that way. Using the conditions 
$^ = 0, $3 = 0, and $4 = the Hamiltonian density can be rewritten as 



9_pa pa I _}_/^a ^a , A"-]"^ _|_ 



a \2 
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where m, n = 1, 2 and the expression contains the physical variables and along 
with the nonphysical A'^. The constraint which is still left relates Aq to the physical 
variables 

5|As+/(z^™p^r = 0. (19) 

Thus, the system is defined by the Hamiltonian density (18) and the constraint (19). 

Let us now turn to the discussion of the spectrum of the system (18-19) which 
is placed in a finite three- volume denoted by Vm- Calculating the spectrum we are 
going to keep only "slow" modes, i.e. the modes with zero momenta but a nonzero 
energy. All the "fast" modes with nonzero momenta can be thought of as being 
integrated out. The net result of the corrections due to the "fast" modes is just a 
perturbative splitting of the energy levels determined by the "slow" modes (for a 
detailed discussion see ref. ^^) Adopting that approximation one can drop all 



the spatial derivatives in the expression for the Hamiltonian and constraint equation 
assuming that all the canonical variables depend on the time variable only. 

Let us turn to the Hamiltonian instead of the Hamiltonian density. Dropping all 
the spatial derivatives one writes down 



d'^^M pa pa , 

2 ^ 4^2 



It is convenient to perform the following rescaling of the canonical variables 

g2/3 1 

The new, rescaled variables are dimensionless. In terms of these variables the expres- 
sion for the Hamiltonian takes the form 

,2/3 



J^YM — — 



Vm 



1/3 



1 1 

_ pa pa \ _ / £abc Ab /ic\2 

2 m m 4^ m-^nj 



(20) 



**The crucial point in this discussion is that the spectrum is calculated in a small volume limit. 
As we mentioned in Section 1, this corresponds to the weak coupling approximation. As a result, 
corrections due to the "fast" modes are of order = {j^^Y^^ @1 and can be neglected in the 
leading approximation. 
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and the constraint equation is given as follows: 



r'"^A^P^ = 0, m = l,2. (21) 

This is the system which defines the spectrum|. The first thing to notice is that 
the potential in eq. (20) has fiat directions. Thus, one would expect a continuous 
spectrum without a mass gap. However, it was proved in ref. (see also p8| ) 



that in the quantum theory, contrary to the naive classical expectation, the operator 
defined in eq. (20) has only discrete positive eigenvalues. As a result, the following 
expression for the spectrum emerges: 

(positive number) 

-f^pcrt = \^YM) = Op^Tt^M = —1/3 ■ 

This expression was used earlier in eq. (17). The exact calculation of the positive 
numbers occurring in the expression above is a complicated problem of YM theory. 
However, as it will be shown below, one can use some analogies and calculate the 
spectrum exphcitly. We turn now to that discussion. 

3.2. The Membrane Matrix Model 



It was shown some time ago ||5T| that the Hamiltonian of a closed bosonic 



membrane in the light-cone gauge can be reduced to the form given in eq. (20). 

The variables substituting the gauge fields in that case occur as coefficients of the 
harmonic expansion of the spatial coordinates on the membrane world surface. The 
two Hamiltonians, one for the membrane and the other one given in eq. (20) formally 
look similar. 

The YM theory constraint (21) also has an analog in the case of the closed mem- 
brane theory. The constraint in that case is related to the residual reparametrization 
invariance of the membrane action which is still left in the light-cone gauge. 

Below we discuss briefiy the membrane action and the way it reduces to the form 
given in eq. (20). Then we deduce the matching condition relating the spectrum of 



§ The operator in (20) acts on functionals of the canonical variable while the momentum operator 



is defined as P„i — "^jjs 
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the closed bosonic membrane to the spectrum of YM theory. The matching condition 
allows one to obtain the spectrum of YM theory by calculating the spectrum of the 
closed bosonic membrane. 

We present below only the basic features of the membrane Hamiltonian construc- 
tion in the light-cone gauge. For details we refer to the original papers [Q, |]51 |. 



The membrane action in flat Minkowski space-time can be written as 

Sm = -T j c/Vy^|det^i,|, (22) 

where T is the membrane tension, the constant with the dimensionality of mass 
cubed; cro,cri,cr2 are the coordinates on the membrane world volume; gij denote the 
components of the induced metric in the membrane world volume 

,..(a) = ^-^^-^ (23) 

where X^, /i = 0, 1, 2, 3, are the space-time coordinates. 

The membrane action is reparametrization invariant. Thus, in accordance with 
the Noether second theorem not all of the variables in the action are independent (as 
in gauge theories). One should carry out the gauge fixing procedure. It is convenient 
to introduce the light-cone coordinates 

X± = i=(X^±X°), 

and choose the light-cone gauge 

X+(a) = X+(0) + (To. 

The light-cone gauge does not completely fix the gauge freedom of the membrane 
action As a result, there still is a residual local invariance left. Hence, one should 
expect to have the Hamiltonian of the theory accompanied by a constraint equation. 
The detailed discussion and the construction of the Hamiltonian is given in refs. ||5( 



51| . We present the final result here. The expressions for the mass squared operator 
and the constraint can be written as follows: 

M2 



2 



1 

_-T)a -pa _|_ /abc-y-b 'y'C\2 

2^' m' m ' m n) 



(24) 



^^As opposed to the case of a string action where in the Hght-cone gauge no freedom is left 

27 



g'^''^XlV'„^ = Q m = l,2. 

The canonical coordinates and momenta are the functions of the time variable only. 
The coordinates in this expression are the coefficients of the harmonic expansion 
of the space-time coordinates on the surface of the membrane. For example, if 
the membrane has the topology of a sphere, then the harmonic expansion mentioned 
above is just the expansion of the space-time coordinates in the basis of spherical 
functions 

oo 

X^{a) = J2 ^2), a= 1, 2...CX), 

a=l 

where y"(cr)'s are the harmonic functions on the sphere. 

If the membrane has the topology of a sphere or torus the harmonic functions 
F"((t) form a representation of the Lie algebra of the SU(oo) gauge group ||51[|p|. 



Thus, the SU{oo) gauge group appears due to the reparametrization invariance of 
the membrane action^. 

The expression (24) resembles the Hamiltonian of the YM system in the approx- 
imation given in eq. (20) and in the 00 limit. The constraint equations in the 
two cases are also similar. 

In order to make use of this analogy let us perform the following rescaling of the 
canonical variables |^ 

p ^ T^l^-p and X T'^/^X. 

The new canonical variables are dimensionless. The expression for the mass squared 
operator in terms of those variables takes the form 



M2 „o/, ri„„„„ 1 



2 



_-T5a -pa I „ibc -y-b vc\2 



(25) 



Thus, one concludes that the spectrum of a closed bosonic spherical membrane is 
determined by the same differential operator as the one for YM theory in the large 
^'''The SU(cx)) group (and its Lie algebra) should be understood as a limit of the SU(N) group at 
TV -> 00. 

§The supersymmetric version of the membrane matrix model is used for the formulation of the 

M theory in the infinite momentum frame JSSt - 

WOne can check that the rescaling procedures we perform lead to dimensionless canonical coor- 
dinates and momenta which satisfy the commutation relation with the unity on its r.h.s. 
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Nc limil||. Matching this expression with eq. (20) one finds the relation between the 
spectrum of YM theory in a finite volume and the spectrum of the closed bosonic 
membrane 



^2/3 /l#2 

- ^1/3 27-2/3' '^^"^ 



where M„'s are the mass eigenvalues defined by the operator given in eq. (25). 
The complimentary constraint equations acting on the physical states ensure that 
the physical eigenfunctionals of the Hamiltonians in both eqs. (20) and (25) are 
the functionals of "colorless" (gauge invariant) variables only. Indeed, in both cases 
the constraint equations (the Gauss's law and its membrane counterpart) serve as the 
generators of the "gauge" transformations of the initial system. Since those generators 
are supposed to annihilate any physical state (imposed as the Gauss's law annihilating 
a state), then all the physical eigenfunctionals should be gauge invariant. 

In the next subsection we calculate for a closed membrane with the topology 
of a sphere and using the matching condition (26) deduce the energy levels for the 
YM theory excitations identifying them with the pseudoscalar glueballs. 

3.3. Calculating the Membrane Spectrum 

We start with a closed spherical membrane. The space-time coordinates on the 
membrane world surface are given as 

= (t, r(t)sm9cos(p, r(t)sm9smip, r{t)cos9), 
<e <TT, < < 27r, 

where r{t) is the time dependent radius of the membrane. 

The induced metric on the membrane worldvolume has the following nonzero 
components: 

gtt = l- r'^it), gee = -r'^(t), g^^ = -r^(t)sin^9. 



^In both cases these differential operators act on the functionals of the canonical coordinates. 
Those functionals are annihilated by the constraint equations. The boundary conditions will be 
fixed later during the actual calculations. 
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The action functional for the membrane takes the form 

S = -T j dtded(pr^{t)sme^l - r^{t). 
Thus, the Lagrangian can be written as follows: 



L{t) = -A'KTr\t)^l-r'^{t). 
Calculating the canonically conjugate momentum 

one derives the Hamiltonian for the spherical membrane^ 

H = + 167r2T2r4. 

As we mentioned above we are looking for the mass squared operator for the mem- 
brane (see eq. (25)). Thus, we need to solve the following Schrodinger equation 

with the boundary conditions ^['(cxd) = and ^E'(O) = 0. 

It is useful to turn to the dimensionless variable z defined as 



z = r 



In terms of z the Schrodinger equation takes the form 



" W)=—r^ + z']^{z). (27) 



(167r2T2)i/3 ^ ^ y dz^ 
The Schrodinger equation with the quartic potential has been extensively studied in 



the literature (for a review see ref. [Q). The results of numerical calculations of the 



first ten eigenvalues can be found in ref. Those calculations are usually done 



for the potential defined on the whole z axis. In our case z is defined on the positive 



W The Hamiltonian looks similar to the one for a relativistic particle with the time-dependent 
mass m(t) — 4TTTr'^{t) and describes the pulsation of the spherical membrane. 
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semiaxis only. Thus, only the odd parity solutions are relevant for the present case. 
Those solutions have nodes at 2; = and satisfy the boundary conditions ^'(O) = 
and ^(00) = 0. 

Here, we present only the first two parity odd eigenvalues of eq. (27) 

° - 2.393644, 



(167r2r2)V3 

Ml 



7.335730. 



(167r2T2)i/3 

Using these expressions and the matching condition (26) one calculates the first two 
energy levels for the YM system in the finite volume 



EU = ^^^^V^2.393644, 



(28) 

o2/3 (167r2)l/3 

E'ym = :,i/3 9 7.335730. (29) 

The numerical values for the energy levels are determined by the strong coupling 
constant g and also by the volume of the domain M.. The strong coupling constant 

is supposed to be taken at the scale appropriate for corresponding glueballs. 

Let us take eqs. (28) and (29) and substitute them into eq. (17). This leads to 

the expression for the total energy inside of the finite volume we are discussing 

1 ?/2 
^EiyM) = f'AxVu + (30) 

where in accordance with eqs. (28-29) 

1^2^^2/3^2712^)1/32393544^ ^^^^^ uj = g"^/^ {27:^^^7.335730. (31) 

The expression (30) can be minimized w.r.t. the value of the three- volume Vm- We 
denote the optimal value for the volume by V, hence ^^^§^f^\v — 0- Using this 
condition and taking the derivative of eq. (30) one finds 

lo'^X^l^,, (32) 



and the value of the total energy for the optimal volume 



AEiVn)=mn--^. (33) 
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Here, we denote by m„ the mass of the corresponding ra'th glueball and by Vn the 
corresponding optimal value of the volume element. Thus, knowing the value of the 
strong coupling constant at the scale appropriate for the lightest glueballs (which is 
about 1.5 — 2.5 GeV) and also knowing the value of the effective size of the YM 0^^ 
glueball state one can predict the value of its mass by means of eq. (33). 

We present below the results of calculations for three different values of the strong 
coupling constant as- The reasonable estimate for the lightest pseudoscalar glueball 
radius is Rq = 0.7 — 1.0 fm ^Of, |25|. The size of the second excited glueball state 
Ri is not known. However, using eqs. (32) and (31) one can estimate that Ri ^ 
1.3i?o = (0.9 — 1.3) fm. The results of numerical calculations of glueball masses for 
those values of the coupling constant and radii are presented below. 
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These predictions can be compared with the result of the lattice calculation for the 



32 



lightest pseudoscalar glueball mass mo = 2.3 ±0.2 GeV [|]. We should stress that the 
masses presented above give just the large approximation to the actual values. We 
regard these numbers as reasonable estimates for the pseudoscalar glueball masses. 

Let us now discuss an interesting consequence of eq. (32). If one knew the effective 
size of the glueball and also the value of Ax, then one would be able (using eq. (32)) 
to calculate the value of the 6 parameter 

In general, the value of Ax is not known. However, in order to get an order of 
magnitude estimate for 9 one can crudely approximate Ax by the lightest glueball 
contribution /gmQ multiplied by the number of glueballs in the spectrum of the 
model (let us call that number A^): Ax ~ Nf^ml ^ A^(200 MeV)^ @, @- Then, 



if as = 0.3 and the lightest glueball radius i?o = 0.8 fm the 6 parameter should be 
equal to ~ 6/ -\/iV- One can also estimate the magnitude of 9 for different values of 
the radius. Generically, if the value of N is not too large, the magnitude of 6 is of 
order of the unity or so. 

Some comments are in order here. First of all the estimate for the 6 parameter 
presented above appears as a result of the physical picture of the glueball formation 
discussed in this work. However, the method of modeling the glueball spectrum by 
means of the membrane Hamiltonian does not depend on a particular mechanism of 
the formation of glueballs. Indeed, whatever the mechanism of the formation is, the 
glueball can always be regarded in some extent as a close domain of space where the 
YM excitations are confined and the spectrum of which is determined by the YM 
Hamiltonian given in eq. (20). 

The second comment concerns the strong CP problem. In this work we deal with 
pure YM theory. No light quark degrees of freedom were included. The large value of 
the 6 parameter that we derived should somehow be neutralized when quark degrees 
of freedom are taken into account. 

There are some possibilities for that. We list below three of them. 

In full QCD the parameter which defines the magnitude of the strong CP violation 
is the sum of the 6 angle used in this work and the phase of the determinant of the 
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quark mass matrix, arg detM. It is possible that those two contributions compensate 
each other and the strong CP violation, being present in pure YM theory, does not 
appear in full QCD. This might lead to an interesting pattern of mixing between the 
pseudoscalar glueball and the rj' meson, when the pseudoscalar glueball being present 
in pure YM model, might not appear in full QCD as a separate state. 

The second possibility is realized if one has a massless quark in the model. In that 
case the 6 dependence can be eliminated from the QCD Lagrangian by an appropriate 
chiral rotation of that quark field. What happens with the glueball state in full QCD 
remains to be studied. 

Finally, one can argue (using the results of refs. |T^, |]17[) that in full QCD the 
rj' meson, mediating interactions between topologically charged objects, provides a 
sufficient (from the experimental point of view) screening of the topological suscep- 
tibility even in the massive theory. In terms of eq. (9) that can be understood by 
including the r]' contribution on the r.h.s. and deducing a Witten-Veneziano type 
relation. More detailed studies of full QCD are needed in order to determine which 
of the above scenarios (if any) can actually be realized. 



Discussions 

In this paper we studied some properties of the YM vacuum which should be 
responsible for the formation of the glueball states. 

The properties of the correlator of the vacuum topological susceptibility as a 
function of the volume element V are discussed. In the weak coupling (small volume) 
approximation it is an increasing function of the argument. Increasing the volume 
continuously the theory passes through a crossover region after which it should be 
regarded as a strongly correlated one. Above the crossover region the topological 
susceptibility becomes a rapidly decreasing function of the argument and reaches its 
asymptotic value (not necessarily zero) in the large volume limit. Thus, the value of 
the vacuum topological susceptibility is screened if the strong coupling regime of the 
theory is considered. 

It is shown that the presence of the 6 angle in the theory along with the screening 
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phenomenon can lead to the formation of a glueball state. An important ingredient 
of that scenario is the existence of the three-form composite field propagating the 
Coulomb-like interaction. 

The spectrum of the YM Hamiltonian resembles in the zero momentum approxi- 
mation the spectrum of a closed bosonic membrane. Using that analogy and calculat- 
ing the spectrum of a closed bosonic membrane we estimate the masses of glueballs 
in the large limit. The result for the lightest glueball is in agreement with 
the lattice prediction. We also predict the mass of the next-to-lightest glueball. This 
result can be checked in future lattice calculations. In general, our approach allows 
us to compute the mass of any heavier glueball state (if such a state exists). The 
method of calculation of the spectrum is in general independent of the mechanism by 
which glueballs are formed in YM theory and the YM vs. the membrane Hamiltonian 
analogy utilized for that calculation can always be applied. 

Notice that the large Nc arguments were not used while deriving eq. (17). The 
large Nc approximation was adopted later on in order to calculate the "positive num- 
ber" occuring on the r.h.s. of eq. (17). Thus, the approach and equations presented 
in this work are not peculiar to the Nc —>■ oo limit. They should rather have some 
wider range of validity beyond the large Nc approximation. For instance, the second 
term in eq. (17) can be thought as a result of the uncertainty principle alone. 

There are a number of interesting questions left out of the discussion in the present 
paper. First of all we did not discuss the fate of a scalar glueball. The effective 
Lagrangian approach to the 0"'""'" channel of pure YM theory was developed in refs. 



41| , . One can apply the YM Hamiltonian vs. the membrane Hamiltonian analogy 
to the calculation of the scalar glueball mass too. This last would correspond to 
the lowest parity-even solution of the Schrodinger equation (27). Hence, the scalar 
glueball would emerge to be lighter than the pseudoscalar one. This is in agreement 
with what is known from various lattice and theoretical studies 0, However, 
the mechanism of the formation of the scalar glueball can not be captured by our 
analysis. 

We did not discuss here how colored degrees of freedom are confined inside of a 
finite closed volume. It was rather assumed that QCD provides this property by some 
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mechanism. Formally, it was assumed that the operator in the Gauss's law, being the 
generator of gauge transformations, should annihilate all the physical states. Thus, all 
those states are supposed to be colorless states by the construction. In various models 
of hadrons, confinement can be warranted by imposing some boundary conditions on 
fields, as in the case of the MIT bag model or the model of ref. or by 



postulating some specific dielectric properties of the vacuum as in the case of the 
Friedberg-Lee model Some discussions of these issues from the point of view of 



QCD can be found in ref. |^ 



Finally, one needs to know what happens when quark degrees of freedom are 
also included in the theory. In that case the mixing between the r/' meson and the 
glueball should play an important role (if those two states exist simultaneously). Our 
discussion of the three-form field in that respect becomes crucial. It is known that the 
r( meson couples to the topological charge density, hence it couples to the three-form 
potential too. Thus, one can naturally couple the 77' meson to the glueball by means 
of the three-form field. These and other related questions will be addressed elsewhere. 



Acknowledgments 

The author is grateful to G.R. Farrar, M. Porrati and M. Schwetz for reading the 
manuscript and for useful suggestions and discussions. I wish to thank T. Banks for 
valuable discussions. The work was supported by grant No. NSF-PHY-94-23002. 

Appendix 

In this appendix we consider the dispersion relation for the correlator of the vac- 
uum topological susceptibility in momentum space. The space-time is assumed to 
be a Euclidean one. The correlator is defined as in eq. (1), Section 1. Before we 
go further let us specify how singularities are handled in eq. (1). The product of 
two operators of the topological charge density is singular as x — * 0. The leading 
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perturbative singularity at x — > can be calculated 



TQ(x)Q(0)oc^ + of^ 



) 



Upon integration in eq. (1) this expression yields a divergent term. A simple way to 
handle the divergence is to allow a small momentum k to flow through the correlator 
function treating x{^) ^ ^lero momentum limit of the corresponding momentum- 
dependent renormalized Green's function 



where k is the regularizing momentum. This relation implies that the limiting pro- 
cedure is supposed to be carried out after the integration and renormalization of the 

divergent parts are already done in momentum space. In what follows we adopt this 
prescription. 

Another type of divergence occurring in eq. (1) is related to the x ^ oo limit. In 
that limit 



Supposing that generically the VEV of the topological charge density might not be 
zero in a CP violating model, one gets the divergence in eq. (1) as y ^ oo. In 
order to eliminate this divergence one can work with the subtracted correlator. This 
amounts to saying that the actual integrand in eq. (1) is the function with the 
following subtraction 



The subtracted function goes to zero in the x ^ oo limit. The coordinate- independent 
subtraction term does not affect our analysis and was dropped for simplicity in Section 
1. It will also be omitted below. 

In what follows we show that continuum contributions vanish in the limit k'^ 0. 

The dispersion relation for the correlator of the vacuum topological susceptibility 
in momentum space can be written as 



/ e'^''{Q\TQ{x)Q{0)\Q)d''x 



-| ren 



X{V) = \ymV^-{k\V) = Ym. 



(o|rg(x)g(o)|o)^(o|g|o)(o|g|o). 



(o|Tg(a;)g(o)|o)-(o|g|o)(o|g|o). 




(Al) 
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where p{s) = Im x(s + ie). The correlator at zero momentum is denoted by x(0)- 
The quantity x'(0) stands for the derivative of the correlator w.r.t. k"^ at k"^ = 0. 

In order to make the integral convergent, and also to account for the correct 
asymptotic behavior of the correlation function at A;^ — ^ oo, we have introduced the 
subtraction terms in the dispersion relation (Al). 

The dispersion relation in the form given is eq. (Al) is not convenient for our 
purposes. In the limit fc^ — * it turns into a trivial identity. One needs to rewrite 
(Al) in a form similar to the one given in eq. (6). For this purpose let us use the 
following identity: 



k^ Ilk 



2 



s^(s — k"^) s — k"^ s 

Substituting this formula into eq. (Al) one rewrites the dispersion relation in the 
following form: 

Go 

where 

1 /■°° p(s)ds , ,, , 1 r°° p(s)ds 



The form of the relation given in eq. (A2) is very formal one. The constants do, bo 
and the integral on the r.h.s. are divergent quantities. When these terms are put 
together all divergences cancel and the whole expression is finite. The divergences 
mentioned above are related to perturbative contributions to the spectral density. 
Thus, it is convenient to separate nonperturbative and perturbative terms. We found 
it useful to apply the decomposition usually adopted in QCD sum rule calculations 
p7| . One decomposes the expression for the spectral density 

p{s)=p'^P{s)+p^\sns-so), 

where the superscripts "np" and "pt" denote nonperturbative and perturbative terms, 
respectively. Here denotes the step function. The constant so sets the continuum 
threshold (or the duality interval) |^| and by the definition so > ""^Gq- assumed 



in this approach that resonance contributions are defined by the nonperturbative 
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part of the spectral density. One also supposes that due to asymptotic freedom 
continuum contributions above the continuum threshold can be approximated by 
leading perturbative terms |^ . 



Let us make the same formal decomposition for the quantities do and 60 

1 pP\s)ds 



do = d+dP\ d^' = / 

71 J so 

bo = b + W\ = - 



^ 1 f<=° ff\s)ds 



vr J so 

Here d and b are the quantities determined by the complicated vacuum structure 
of YM theory. As we mentioned already, in the weak coupling approximation with 
noninteracting instantons d is defined as the value of the topological susceptibility 
of a dilute instanton gas in the large volume limit of pure YM theory. The quantity 
d appears in eqs. (5) and (6-8) in the text. Finally, using all the expressions given 
above one derives 

We should notice here that eqs. (A3) and (Al) differ from each other by some formal 
redefinitions. Moreover, eq. (A3) is written adopting some particular scheme of 
separation between perturbative and nonperturbative contributions. That procedure 
is not unambiguous. In that respect, eq. (A3) should be regarded as an expression 
defined within the framework of the particular prescription outlined above. 

Now one can use the fact that the quantity k"^ is a regularizing momentum. Thus, 
one can assume that k"^ is very small, so that the condition sq >> k"^ is readily 
satisfied. The last integral on the r.h.s of eq. (A3) can be expanded in a power series 
of the ratio fc^/so (since that integral is convergent). Performing the expansion, and 
then Fourier transforming eq. (A3) with the weight -^2:^; one derives the expression 
for the correlator {Q\TQ{x)Q{Q)\Q) in the following formp] 

(0|TQ(a;)Q(0)|0) = d 6^^\x) - b d^6^^\x) + - / p"P(s)Di.(v^|s|) ds+ 

Go 

^We use the following normalization for the delta function: = j^py^ ^^^^ df^k and 
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vr s \s J 

Eq. (A4) is a general form of the expression given in eq. (6) in Section 1. In order to 
reproduce the sum on the r.h.s. of eq. (6) one needs to make the following substitution 
in eq. (A4) 

n 

Some terms on the r.h.s. of eq. (A4) with derivatives of the Dirac delta function yield 
vanishing contributions upon integration in eq. (1). For that reason those derivative 
containing terms were omitted in eq. (6). 
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